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Abstract 

The approximate analytic bound state solutions of the Klein-Gordon equation with equal scalar 
and vector exponential-type potentials including the centrifugal potential term are obtained for 
any arbitrary orbital angular momentum number / and dimensional space D. The relativistic/non- 
relativistic energy spectrum equation and the corresponding unnormalized radial wave functions, 
in terms of the Jacobi polynomials P^'^\z), where a > —1, /? > —1 and z G [— 1,+1] or the 
generalized hypergeometric functions 2-^1(0^) b; c; z), are found. The Nikiforov-Uvarov (NU) method 
is used in the solution. The solutions of the Eckart, Rosen-Morse, Hulthen and Woods-Saxon 
potential models can be easily obtained from these solutions. Our results are identical with those 
ones appearing in the literature. Finally, under the PT-symmetry, we can easily obtain the bound 
state solutions of the trigonometric Rosen-Morse potential. 
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I. INTRODUCTION 



The exact solutions of the wave equations (non-relativistic or relativistic) are very impor- 
tant since they contain all the necessary information regarding the quantum system under 
consideration. However, analytical solutions are possible only in a few simple cases such 
as the hydrogen atom and the harmonic oscillator [1,2]. Most quantum systems could be 
solved only by using approximation schemes like rotating Morse potential via Pekeris ap- 
proximation [3] and the generalized Morse potential by means of an improved approximation 
scheme [4]. Recently, the study of exponential- type potentials have attracted much atten- 
tion from many authors [5-26]. These potentials include the Woods-Saxon [5,6], Hulthen 
[7-16], Manning-Rosen [17-22], the Eckart [23-25] and the Rosen- Morse [26] potentials. 

The spherically symmetric Eckart-typc potential model [27] is a molecular potential model 
which has been widely applied in physics [28] and chemical physics [29,30] and is generally 
expressed as 

V{r; q) = Vi cos ec/i^ar — V2 coth^ ar, Vi, V2 > 0, — 1 < g < or g > 0, (1) 

where the coupling parameters Vi and V2 describe the depth of the potential well, while 
the screening parameter a is related to the range of the potential. It is a special case of 
the five-parameter exponential- type potential model [31,32]. The range of parameter q was 
taken as 5 > by Ref. [33] and has been extended to— l<g<0org>0or even complex 
by Ref. [34]. The deformed hyperbolic functions given in (1) have been introduced for the 
first time by Aral [35] for real q values. When q is complex, the functions in (1) are called 
the generalized deformed hyperbolic functions. The Eckart- type potentials (1) can also be 
written in the exponential form as 

-2ar I I ^g-2ar 

Vir] q) = 4Vi ^ - V2 (2) 

The study of the bound and scattering states for the Eckart-type potential has raised a great 
deal of interest in the non-relativistic as well as in relativistic quantum mechanics. The s- 
wave (/ = 0) bound-state solution of the Schrodingcr equation for the Eckart potential has 
been widely investigated by using various methods, such as the supersymmetric (SUSY) 
shape invariance technology [36], point cannonical transformation (PCT) method [37] and 
SUSY Wentzel- Kramers- Brillouin (WKB) approximation approach [38]. The bound state 
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solutions of the s-wave Klein-Gordon (KG) equation with equally mixed Rosen-Morse-type 
(Eckart and Rosen-Morse well) potentials have been studied [39] . The bound state solutions 
of the s-wave Dirac equation with equal vector and scalar Eckart-type potentials in terms of 
the basic concepts of the shape- invar iancc approach in the SUSYQM have also been studied 
[24] . The spin symmetry and pseudospin symmetry in the relativistic Eckart potential have 
been investigated by solving the Dirac equation for mixed potentials [25]. Unfortunately, 
the wave equations for the Eckart-type potential can only be solved analytically for zero 
angular momentum states because of the centrifugal potential term. Some authors [23- 
25] studied the analytical approximations to the bound state solutions of the Schrodinger 
equation with Eckart potential by using the usual existing approximation scheme proposed 
by Greene and Aldrich [40] for the centrifugal potential term. This approximation has 
also been used to study analytically the arbitrary /-wave scattering state solutions of the 
Schrodinger equation for the Eckart potential [41,42]. The same approximation scheme for 
the spin-orbit coupling term has been used to study the spin symmetry and pseudospin 
symmetry analytical solutions of the Dirac equation with the Eckart potential using the 
AIM [43]. Overmore, the pseudospin symmetry analytical solutions of the Dirac equation 
for the Eckart potential have been found by using the SUSY WKB formalism [44]. Very 
recently, for the first time, the approximation scheme for the centrifugal potential term has 
also been used in [45] to obtain the approximate analytical solution of the KG equation for 
equal scalar and vector Eckart potentials for arbitrary Z-states by means of the functional 
analysis method. 

This approximation for the centrifugal potential term [7,14,40] has also been used to solve 
the Schrodinger equation [7,14], KG [8,15] and Dirac equation [15] for the Hulthen potential. 
Recently, the KG and Dirac equations have been solved in the presence of the Hulthen 
potential, where the energy spectrum and the scattering wave functions were obtained for 
spin-0 and spin-| particles, using a more general approximation scheme for the centrifugal 
potential [15]. They found that the good approximation, however, occurs when the screening 
parameter a and the dimensionless parameter 7 are taken as a = 0.1 and 7=1, respectively, 
which is simply the case of the usual approximation [7,14]. Also, other authors have recently 
proposed an alternative approximation scheme for the centrifugal potential to solve the 
Schrodinger equation for the Hulthen potential [46]. Taking ou — 1, their approximation can 
be reduced to the usual approximation [7,14]. Very recently, we have also proposed a new 



3 



approximation scheme for the centrifugal term [9] . 

The Nikiforov-Uvarov (NU) method [47] and other methods have also been used to solve 
the £)-dimensional Schrodinger equation [48] and relativistic £)-dimensional KG equation 
[49], Dirac equation [4,10,26,50] and spinlcss Salpctcr equation [51]. 

The aim of this work is to employ the usual approximation scheme [40,45] in order to 
solve the D-dimensional radial KG equation for any orbital angular momentum number I 
for the scalar and vector Eckart-type potentials using a general mathematical model of the 
NU method. This offers a simple, accurate and efficient scheme for the exponential-type 
potential models in quantum mechanics. We consider the following relationship between the 
scalar and vector potentials: V{r) = Vq + PS(r), where Vq and /3 are arbitrary constants 
[52]. Under the restriction of equally mixed potentials S{r) = V{r), the KG equation turns 
into a Schrodinger-like equation and thus the bound state solutions are very easily obtained 
through the well-known methods developed in the non-relativistic quantum mechanics. It is 
interesting to note that, this restriction include the case where V{r) — when both constants 
vanish, the situation where the potentials are equal {Vq — 0; P — 1) and also the case where 
the potentials are proportional [53] when Vq = and /3 = ±1, which provide the equally- 
mixed scalar and vector potential case V{r) = ±S{r). Very recently, we have obtained an 
approximate analytic solution of the KG equation in the presence of equal scalar and vector 
generalized deformed hyperbolic potential functions by means of parameteric generalization 
of the NU method. Furthermore, for the equally-mixed scalar and vector potential case 
V{r) — ±S{r), we have obtained the approximate bound state rotational- vibrational (ro- 
vibrational) energy levels and the corresponding normalized wave functions expressed in 
terms of the Jacobi polynomial Pn^'''\x), where ^ > —1, i> > —1 and x e [—1,-1-1] for a 
spin-zero particle in a closed form [54]. 

The paper is structured as follows: In section 2, we derive a general model of the NU 
method valid for any central or non-central potential. In section 3, the approximate analyt- 
ical solutions of the D-dimensional radial KG equation with arbitrary Z-states for equally- 
mixed scalar and vector Eckart-type potentials and other typical potentials are obtained by 
means of the NU method. Also, the exact s-wave KG equation has also been solved for 
the Rosen-Morse- type potentials and other typical potentials. The relative convenience of 
the Eckart-type potential (Rosen- Morse- type potential) with the Hulthen potential (Woods- 
Saxon potential) has been studied, respectively. We make some remarks on the energy 
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equations and the corresponding wavefunctions for the Eckart and Rosen-Morse well poten- 
tials in various dimensions and their non-relativistic limits in section 4. Section 5 contains 
the summary and conclusions. 



II. NU METHOD 

The NU method is briefly outlined here and the details can be found in [47]. This method 
was proposed to solve the second-order differential wave equation of the hypergeometric-type: 

a\z)ij':,{z) + a{z)T{z)ij'^{z) + a{z)ij^{z) = 0, (3) 

where a{z) and d^{z) are at most second-degree polynomials and t{z) is a first-degree poly- 
nomial. The prime denotes the differentiation with respect to To find a particular solution 
of Eq. (3), one can decompose the wave function ipni^) as follows: 

fPniz) ^(f>n{^)yn{z), (4) 

leading to a hypergeometric type equation 

a{z)y';{z) + T{z)y'^{z) + Xy^iz) = 0, (5) 

where 

X^k + 7r'{z), (6) 
and yn{z) satisfies the Rodrigues relation 

A rl^ 

In the above equation, is a constant related to the normalization and p{z) is the weight 
function satisfying the condition 

a{z)p'{z) + {a'{z)-T{z))p{z) = 0, (8) 

with 

T{z)=T{z)+27r{z),T'{z) <0. (9) 

Since p{z) > and a{z) > 0, the derivative of t{z) should be negative [47] which is the 
essential condition for a proper choice of solution. The other part of the wavefunction in 
Eq. (4) is defined as 

a{z)cf>'{z)-7T{z)ci>{z)^0, (10) 



where 

= I K(^) - r{z)] ± [a'{z)-T{z)f-a{z)+ka{z). (11) 

The determination of k is the essential point in the calculation of 7t{z), for which the dis- 
criminant of the square root in the last equation is set to zero. This results in the polynomial 
Tr{z) which is dependent on the transformation function z{r). Also, the parameter A defined 
in Eq. (6) takes the following form 

X^^ -nr'iz) -]-n{n-l)a"{z), n = 0,l,2,---. (12) 

We may construct a general recipe of the NU method valid for any central and non-central 
potential. We begin by comparing the following hypergeometric equation 

[z (1 - c^z)f i,"Sz) + (1 - c^z) (ci - c^z)\ ij'^iz) + {-Az^ + BZ-C) ^ J^) - 0, (13) 

with its counterpart Eq. (3), we then obtain [54] 

r{z) = ci - C2Z, a{z) = z {1 - c^z) , a{z) = -Az^ + Bz-C. (14) 

Substituting Eq. (14) into Eq. (11), we find 

7r{z) = C4 C5Z ± [(c6 - C3A;+ _) + {c-j + A;+ _) z + , (15) 

where 

1 1 

C4 = - (1 - ci) , C5 = - (c2 - 2C3) , Cq^cI + A, C7 2C4C5 - B, cg = C4 C. (16) 

The discriminant under the square root sign must be set to zero and the resulting equation 
must be solved for /c, it yields 

A;+ _ = - (c7 -h 2C3C8) ± 2^0^, (17) 

where 

cg = C3 {c^ + C3C8) + ce. (18) 

Inserting Eq. (17) into Eq. (15) and solving the resulting equation, we make the following 
choice of parameters: 

7r(2;) = C4 c^z - [{sfc^ + c^y/c^) z - ^/c^ , (19) 
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A;_ = - (C7 + 2c3C8) - 2y^. (20) 

Further, from Eq. (9), we get 

r{z) = 1 - (c2 - 2c5) z-2 [(^ + C3x/^) z - , (21) 

whose derivative must be negative: 

r'{z) = -2c3 - 2 (^/^ + caV^) < 0, (22) 

in accordance with essential requirement of the method [47]. Solving Eqs. (6) and (12), we 
get the energy equation: 

(C2 - ca) n + can^ - (2n + 1) Cg + (2n + 1) ( + ca ^c^) + C7 + 2C3C8 + 2^/^ = 0, (23) 

for the potential under investigation. Let us now turn to the wave functions. The solution 
of the differential equation (8) for the weight function p{z) is 

p{z)=z''°{l-cszy'\ (24) 

and consequently from Eq. (7), the first part of the wave function becomes 

yn{z) = Pi''"''''\l-2c3z), Re(cio)>-l, Re(cii) > -1, (25) 

where 

2 

ClO = Ci + 2C4 + 2^/c^ - 1, Cii = 1 - Ci - 2c4 H \/c^, (26) 

C3 

and Pi"'''^(l — C32;) are Jacobi polynomials. The second part of the wave function (4) can be 
found from the solution of the differential equation (10) as 

0(^) = ^'^l^(l-C3^)^^^ (27) 

where 

C12 = C4 + v^, ci3 = -C4 + — (v^ - C5) . (28) 

C3 

Hence, the general wave functions (4) read as 

ui(z) = Nr^z'^^'il - C3^)'^"Pi'^^°'^")(l - 2C3^), (29) 
where is a normalization constant. 



III. BOUND-STATE SOLUTIONS 



The £)-dimensional time-independent arbitrary Z-states radial KG equation with scalar 
and vector potentials S{r) and V{r), respectively, where r = |r| describing a spinless particle 
takes the general form [3,49]: 

VL^£.7^5(x) + ^ { [K. - V{r)r - [Mc' + S{r)Y} i^tlZ'^^) = 0' 

= E ^' W^^r;^ W = Riir)Yi^!,,je„ e,,.--, On-i), (30) 

j=l 3 

where Eni, M and denote the KG energy, the mass and the D-dimensional Laplacian, 
respectively. In addition, x is a D-dimensional position vector. Let us decompose the radial 
wave function Ri{r) as follows: 

i?,(r) = r-(^-^)/2^,(r), (31) 

we, then, reduce Eq. (30) into the £)-dimensional radial Schrodinger-hke equation with 
arbitrary orbital angular momentum number / as 

^ + {l^-- - ''(^)l^ - [^^^ + S(r)f - M!£!} „,„ . 0, (32) 

where we have set = [{M - 2f - 1] /4 and M = D+21 where ^ = 0, 1, 2, • • • . Under 

the equally mixed potentials S{r) = ±V{r), the KG turns into a Schrodinger-like equation 
and thus the bound state solutions are very easily obtained with the help of the well- 
known methods developed in the non-relativistic quantum mechanics. We use the existing 
approximation for the centrifugal potential term in the non-relativistic model [7,14] which 
is vahd only for g = 1 value [49,55]: 



V{r) = ^^i^ « Aa^ni' + 1)^^^^^^, I'^iM- 3) /2, (33) 



in the limit of small a and I'. 
A. The Eckart-type model 

At first, let us rewrite Eq. (2) in a form to include the Hulthen potential. 
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V{r; q) = W-, ^ - V2- ^ - V^-^ —, (34) 

and then follow the model used in Refs. [49,55,56] by inserting the above equation and the 
approximate potential term (33) into (32), we obtain 

+ r. — } Mr) 



h^c^ \ (1 - qe-'^^^f (1 - ge-2"^) 

^ ^(Mc2)'-£;2jxi,(r), kKO) = 0, (35) 



which is now amenable to the NU solution. We further use the following ansatze in order to 
make the above differential equation more compact 



z{r) = e-^-, = ^ /3 = ^ ' ' + I' {I' + 1), 

7 = , A = — , Q = 2hca. (36) 

Notice that \Eni\ < Mc^. The KG equation can then be reduced to 



[z{l - qz)Y + z{l - qzY 



dz"^ dz 
+ {-q'iel, + X)z' + {2qel, + qX - qj - P)z - {el, - 7) } Ui{z) = 0, (37) 

where r G [0, 00) — )■ 2; G [0, 1]. Before proceeding, the boundary conditions on the radial 
wave function ui{r) demand that ii;(r — )■ 00 or 2; ^ 0) ^ and M/(r = or 2; = 1) is finite. 
Comparing Eq. (37) with Eq. (13), we obtain values for the set of parameters given in 
section 2: 

Ci = 1, C2 = C3 = g, C4 = 0, C5 = -|, C6 = q^ (^li + ^ + 1^ ' 



C7 = -g ( 24 + A - 7 - -j , C8 = - 7, C9 = (2 j ^ V 1 ' " ~ ^' 



4/9 _ z:^ — : _ 1 4/9 



cii = ^/ 1 + y , C12 = ^^el^ - 7, ci3 = - ^^1 + y 1 + ^ 
^ = ?'(4 + A), 5 = g(24 + A-7-^), ^^ = 4-7, (38) 



and the energy equation via Eq. (23) as 



.il+Af (n + 5£ 7^ n-0 12... (39) 
where 5 — \ y^l + . Making use of Eq. (36), the above equation turns to become 

M^c^ - El, = {hcaf in + wf + ^ "^f ^ .^^ + ^f,' + ± Mc^) (l^, - V^) , (40) 

(zftca) [n + w) 

where w = ^ (^1 + ,^Jl+ ^^'^^^ + ~ ^^ll^,caf ^ "^) ' "^^^ energy Eni is defined imphcitly by Eq. 
(40) which is a rather comphcated transcendental equation having many solutions for given 
values of n and I. In the above equation, let us remark that it is not difficult to conclude 
that bound-states appear in four energy solutions; only two energy solutions are valid for 
the particle = E^, and the second one corresponds to the anti-particle energy E'^ = E~i 
in the Eckart-type field. 

Referring to the general parametric model in section 2, we can also calculate the corre- 
sponding wave functions. The explicit form of the weight function becomes 



p{z) = z'^{l-qzr--\ p=\ 



{Enl±M^){V2 + V^) 1 

n + w - — 

2 {hca) n + w 



(41) 



which gives the following Jacobi polynomials: 

yn{z)^Pi'^''--'\l-2qz), (42) 

as a first part of the wave functions. The second part of the wave functions can be found as 

0(z) ^ z^il - qzy. (43) 

Hence, the unnormalized wave functions expressed in terms of the Jacobi polynomials read 

ui{z) = Unz^{l - g^)'"Pi^f'2"'-^)(l - 2?^), (44) 

and consequently the total radial part of the wave functions expressed in terms of the 
hypergeometric functions are 

Ri^r) = A/;r-(-°-^)/2 (e-2"^7 (1 - ge-^"'')'" ^Fi {-n, n + 2 (p + w) ; 2p + 1; ge-^""), (45) 

where Afn is a constant related to the normalization. The relationship between the Jacobi 

polynomials and the hypergeometric functions is given by Pi"'^'*(l — 2qx) =2 Ei{—n, n -|- a -|- 

00 

6 + 1; a + 1; X), where .F.ii., 7; x) = J] ^^^^^gSf^tf ■ 



fe=0 
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Now, when taking V2 — Vs, the energy equation (40) satisfying Eni for the equally- mixed 
scalar and vector Eckart-type potentials becomes 

M^c^ - El, = {n.af {n + wf + ^ "f^'^' f (46) 

[hca) [n + w) 



and the wave functions: 
ui{z) = Xnz'^il - qz)^'Pi'''^''^-'\l - 2qz), v = \ 



{E^i±Mc')V2 1 
n + w — 



(hca) n + w^ 



(47) 



or the total radial wave functions in (30) are 

Ri(r) = A/;r-(^-^)/2 ^e-2ary _ qe-2ar^w n + 2 {v + w) ■2v + l' ge-""'-), (48) 

where Mn is a normalization factor. The results given in Eqs. (46) and (47) are consistent 
with those given in Eqs. (15) and (18) of Ref. [45]. 

Also, in taking q — 1, 2a ^ a, Vi — V2 — Q and Vz — Vb, Eq. (34) turns to become the 
Hulthen potential. Hence, we find bound state solutions for equally-mixed scalar and vector 
S{r) — V{r) Hulthen potentials in the KG theory with any orbital angular momentum 
quantum number I and an arbitrary dimension D, 

/^^;^^(fa)(« + .)_(Mc^ + g„,)ro 1 D + 21-1 

V 2 hca {n + vY 2 ' ^ ' 

u,iz) = A4 (e-)^ (1 - e-TPi^-^^-Hl " 2z), . = ^ - iM^^Ji^^^ (50) 

^ [hca) n + V 

and the Jacobi polynomial in the above equation can be expressed in terms of the hyperge- 
ometric function: 

Ri{r) = A4r-(^-^)/' (e-"'-)' (1 - e""^)'^ ^F^ (-n, n + 2 (<j + i/) ; 2^ + 1; e""''), (51) 

where Mn is a constant related to the normahzation. The above results are identical to those 
found recently by Refs. [49,57]. 

In the non-relativistic limit, inserting the equally mixed Eckart-type potentials (1) into 
the Schrodinger equation gives 

(fui{r) , / 2ME^i _ W-h'-J'jr -rl)\t-'-"''' 2MVAl + qi:-^n \ . N n 

dr^ ^\ h^[l- qe-^-rf + h^ [l - qe-^ar^ j '''^'"^ " ^' ^^^^ 

and further making use of the following definitions: 

,, = ^E^. B„,<o,;, = ?^ + r(r + i),, = H^.T = 2te, (53) 
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lead us to obtain the set of parameters and energy equation given before in Eqs. (38) and 
(39) with 7 = A. Incorporating the above equation and using Eq. (39), we find the following 
energy eigenvalues: 



nl 



2M 



2t/2 



.2 2/^ X2^M^^ 

n a [n + Wi) + 



fi^a"^ {n + wif 



(54) 



In addition, following the procedures indicated in Eqs. (41)-(45), we obtain expressions for 
the radial wave functions: 



2ar\ 



n + Wi — 



MV2 1 
h?a'^ n + wi 



(55) 



B. The Rosen- Morse- type model 

Under the replacement of q by —q, the Eckart-type potential model given in Eq. (1) will 
become the Rosen-Morse-type potential model given in Eq. (2) of Ref . [39] : 



V{r, q) = Vi sec h^ar — V2 tanhg ar, Vi, V2 > 0, 



or alternatively [26,58] 



V{r, q) = AVi- 



,—2ar 



-2Qr^2 



1 — qe 



—2ar 



We may rewrite the above equation in a form to include the Woods-Saxon potential, 

„— 2ar 1 „„— 2or 



(56) 



(57) 



V{r, q) = 4Vi 



(1 + qe-'^'^^y 1 + qe 



1 + qe 



-2ar ' 



Using the following definitions 



Q 



2(E„,o±Mc2)T/3 



7 = ^ 0) = '—^^ , A - A(Z ^ 0) ^2 

we write the s-wave KG equation with S{r) — ±V{r) for the potential (58) as 



(58) 



(59) 
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+ I (<o + x)z' + qi^-X- 24,0 - f j ^ - «o - 7) I «n(^) = 0. (60) 

Following same procedures used in the previous subsection, we obtain values for the param- 
eters given in section 2: 



ci = 1, C2 = C3 = -g, C4 = 0, C5 = -, C6 = ? ( £„,o + + 4 ) ' 



4/9 



Cll = --^Z 1 - y , C12 = \/ 4,0 - 7, Ci3 = 5 = - I 1 - W 1 - 



~ 1 

2 
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^ = (4,0 + a) , = -g ( 2el, + ^ + A- 7),C = £^ „ _ 



(61) 



and the energy equation 



=-2 _ 



(7 + a)' 



+ 



n + S 



+ 



7 — A 



4(n + 5)2 4 2 

Inserting Eq. (59) in the above equation, we obtain energy equation satisfying E^fi, 



(62) 



2 , (^„,o ± Mc2)' (\/2 + Vs 



{2hcaf' (n + wf' 



1 / _ / 8(K,o±Mc2)\/i 



(63) 



2 \ \/ q (hcaY 

The corresponding unnormalized wave functions can be calculated as before, the explicit 
form of the weight function becomes 



p{z) ^ - qzr^-\ 
which gives the Jacobi polynomials: 



n-\-w 



2 {hcaf n + w 



(64) 



yn{z)^Pi'^''''-'\l + 2qz), 



(65) 



as the first part of the wave function. The second part of the wave function can be found as 



(f){z) zP{l + qzf 



(66) 
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The unnormalized wave function reads 

Un{z) = XnZ^il + qzfPi^^^^''-'\l + 2qz), (67) 

and thus the total radial part of the radial wave functions in (30) can be expressed in terms 
of the hypergeometric functions as 

Rn{r) = Nn (e-2"-)^(l + ge-2"'-)- (-n, n + 2 (p + w) ; 2p + 1; -ge'^"'^), (68) 

where Mn is a normalization factor. 

In taking V2 = V3 in Eq. (63), we find the equation for the potential in (56) satisfying 
Enfl in the s-wave KG theory, 

Mc -E o={hca) [n + w) + ——^ — — 372> (69) 

[nca) [n + w) 

and the wave functions are 

Un{r) = Mn (e-^"")"^' (1 + ge-2«')'2p^2pi,2i5-i)^^ ^ 2qe-^"''), 
= Nn (1 + qe-^^'-f ^F, (-n,n + 2 {p, + w) ; 2p, + 1; -ge'^"'^). 



1 



, „ (K,0±Mc2)l^2 1 

n + w — 



9 ~ , (70) 

{hcaf n + w_ ' 

where A'^ is a normalization constant. After making appropriate change of the potential 
parameter Vi — > —Vi in Eq. (56), our results in Eqs. (69) and (70) become identical with 
Eqs. (13) and (14) of Ref. [39]. 

Also, taking q = 1, 2a ^ a, Vi = V2 = and = —Vq, Eq. (58) turns to become the 
Woods-Saxon potential. Hence, we can find bound state solutions in the s-wave KG theory 
with equally- mixed scalar and vector S{r) = V{r) for Woods-Saxon potentials as 

n , (Mc^ + Ko)Vbl 



M^c^ - El, = hcap2, P2 = -+ ,r (71) 

^ [nca) n 

and wave functions: 

Un{r) = Nn {e-'^T ^F'-'Hl + 2e-"^), (72) 
or alternatively, it can be expressed in terms of the hypergeometric function as 

Rn{r) = Nnv-' {e-"'Y' 2F1 i-n, n + 2p2; 2p2 + 1; e""''), (73) 
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where Mn is a constant related to the normahzation. Under appropriate parameter replace- 
ments, we obtain the non-relativistic limit of the energy eigenvalues and eigenfunctions of 
the above two equations are 



En,0 - 



2M 
and 



nha 2MVq1 

\ ^- 

2 ha n 



, n ^ 0, (74) 

Un{r)^Nn{e--T 2Fi{-n,n + 2p,;2p2 + l;e-n, Pa = ^ + (75) 

^ [nca) n 

respectively, which is the solution of the Schrodinger equation for the potential E(r) = 
V{r) + S{r) = 2l^(r).The above results are identical to those found before by Ref. [6]. 

IV. DISCUSSIONS 

In this section, at first, we choose appropriate parameters in the Eckart-type potential 
model to construct the Eckart potential, Rosen-Morse well and their PT-symmetric versions, 
and then discuss their energy equations in the framework of KG theory with equally mixed 
potentials. 

A. Eckart potential 

Taking g = 1, the potential (1) turns to the standard Eckart potential [27] 

V{r) = Vi cos ech?ar - V2 coth ar, Vl, ^2 > 0. (76) 

In natural units {h — c — l),we can obtain the energy equation (46) for the Eckart potential 
in the three-dimensional spinless KG theory as 



M' - = a'{n + w'Y + 



a' {n + w'Y 

which is identical with those given in Eq. (22) of Ref. [39] under the equally-mixed potential 
restriction given by S{r) = ±V{r). The unnormalized wave function corresponding to the 
energy levels is 

Ri{r) = AC/r-^^-^)/' {e-^'^'-y (1 - e-2-)-'p^2.,2«,'-i)(^ _ ^^-2ar^^ (73) 
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where M^i is a normalization factor. 



(i) For s-wave case, the centrifugal term ^-z^"*"^^ ~ ^ consequently the ap- 

0, too. Thus, the energy eigenvalues 



proximation term {D + 2l-l){D + 2l-i)a^ ^ 
take the following simple form 



, Wi 



2'^ + 



8(K,o±M)l/i 



1 + 



(79) 



(n + wiY 

(ii) In the non-relativistic approximation of the KG energy equation (potential energies small 

compared to Mc^ and E ^ Mc^) Eq. (32) reduces into the form [59] 



2M dr^ 



ui(r) = {E-Mc'')ui(r). 



(80) 



When V{r) — S{r), the energy eigenvalues obtained from Eq. (80) reduces to those energy 
eigenvalues obtained from the solution of the Schrodinger equation for the sum potential 
E(r) = 2V{r). In other words, the non-relativistic limit is the Schrodinger-like equation 
for the potential 8^1-^-^-7 - 21/27^7^. This can be achieved by making the parameter 

(l-e ) J- e 

replacements M + Eji —f 2M and E^ — M E^r, so the non-relativistic limit of our results 
in Eq. (46) reduces to 

1 



E 



NR 



2M 



o? in + + 



2M^Vi 
o? in -|- w-if' , 



(81) 



and the corresponding wave functions in (48) become 

i?;(r) = (1 - e-2°-)"'2p^2t.2,2«,2-i)^i _ 2e-2«'-), 

2My2 1 



V2 = 



n W2 - 



a.^ n -|- W2 



1 
2 



1 + W(1 + 2Z')' + 



a" 



(82) 



The above two equations arc identical with the NU solution of the Schrodinger equation for 
a potential V{r^ (cf. Eqs. (54) and (55)). 



B. PT-symmetric Trigonometric Rosen-Morse (tRM) potential 

When we make the transformations of parameters as a — > ia, V2 — > —iVi-, and V\ — > — T^, 
and using the relation between the trigonometric and the hyperbolic functions sin [iax) — 
isinh [ax] , the potential (1) turns to become the PT-symmetric tRM potential [60]: 

TT 

V{x) = Vx csc^ ax - V2 cot ax, Re{Vi) > 0, a = — , a; = (0, d], (83) 

Zidi 
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where Vi — a{a + 1) and V2 — 2b. This potential is displayed in Figure 1 which is nearly 
linear in 7r/3 < ax < 27r/3, Coulombic in 7r/90 < ax < tt/SO and infinite walls at and tt. 
So it might be a prime candidate for an effective QCD potential. For a potential V{x), when 
one makes the transformation of a; — )> — x and i — )■ —i, if the relation V{—x) = V*{x) exists, 
the potential V{x) is said to be PT-symmetric, where P denotes parity operator (space 
reflection) and T denotes time reversal [6,61]. Our point here is that V{x) interpolates 
between the Coulomb-and the infinite wall potential [62] going through an intermediary 
region of linear-x-and harmonic-oscillator x^ dependences. To see this it is quite instructive 
to expand the potential in a Taylor series which for appropriately small x takes the form of 
a Coulomb-like potential with a centrifugal-barrier like term, provided by the csc^ ax part 
[63], 

V{x) ^ -— + ax < 1. (84) 

Oix [ax) 

For ax < Tr/2 we can then take the potential (84) plus a linear like perturbation 

AV{x) = Vi/3 + V2X/3, (85) 

as an approximation of tRM potential. The potential (83) obviously evolves to an infinite 
wall as ax approaches the limits of the definition interval < ax < n, due to the behavior 
of the cot ax and esc ax for Vi > 0. The potential is essential for the QCD quark- gluon 
dynamics where the one gluon exchange gives rise to an effective Coulomb-like potential, 
while the self gluon interactions produce a hnear potential as established by lattice QCD 
calculations of hadron properties (Cornell potential) [64] . Finally, the infinite wall piece of 
the tRM potential provides the regime suited for the asymptotical freedom of the quarks. 
Now, making the corresponding parameter replacements in Eq. (46), we end up with real 
energy equation for the above PT-symmetric version of the Eckart-type potential^ in the 
KG equation with equally mixed potentials, 

{Mc^y - El, = ^ 5 . , - {hcaf in + wf , (86) 

[nca) [n + w) 

and the radial wave functions build up as 

Rl{x) = MnX-^""-'^^^ ^^+^2axy _ ^+^2a.y n + 2 {v + w) ; 2v + 1; e+*'""), 

1 



^=2 



{Eni±M^)V2 1 

n + w + I 5 

[hca] n + w 



(87) 
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C. Standard Rosen-Morse well 



Taking q^l,Vi^ -Vi {(3 -(3) and V2 -V2 (7 -7), the potential (56) turns to 
the standard Rosen-Morse well [26,58] 



V{r) = -Vi sec h'^ar + V2 tanh ar, Vi, V2 > 0. 



(88) 



This potential is useful in discussing polyatomic molecular vibrational energies. An example 
of its apphcation to the vibrational states of NH^ was given by Rosen and Morse in [26,58]. 
Making the corresponding parameter replacements in Eq. (69), we obtain the energy equa- 
tion for the Rosen-Morse well in the s-wave KG theory with equally mixed potentials. 



- El, 



(89) 



The unnormalized wave function corresponding to the energy levels is 

,(r) = (e-^-)^^ (1 + e-2-)'^^i''"'''"'V + 2e-'"0, 

(K,o ±M)V2 1 



Ur 



Vl = 



n + Si + 



(90) 



a n + Si. 

where Af^ is a normalization constant. The results given in Eqs. (89) and (90) are consistent 
with those given in Eqs. (19) and (20) of Ref. [39], respectively. The s-wave energy states 
of the KG equation for the Rosen-Morse potential are calculated for a set of selected values 
parameters in Table 1. 

When Vo = Sq, the non-relativistic hmit is the solution of the Schrodinger equation for 
the potential —SVi ^^^^l^lrf + l^t-lll ■ In the non-relativistic limits, the energy spectrum 



IS 



Enr — 



2M 

and the wave functions are 



a 



(n + 62^ 



■l+e- 



2 4M2\/2 
+ ^ 



a' 



(n + ^2) 



,S2^-^{1 



1 + 



I6MV1 



(91) 



n + 62 + 



2MV2 1 

n + S2. 
(92) 
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V. CONCLUSIONS 



We have used a parametric generalization model derived from the NU to obtain the ana- 
lytic bound state solutions of the KG equation with any orbital angular momentum quantum 
number I for equally mixed scalar and vector Eckart-type potentials. These calculations in- 
clude energy equation and the unnormalized wave functions being expressed in terms of 
the Jacobi polynomials or the hypcrgeometric functions. Furthermore, making appropriate 
changes in the Eckart-type potential parameters, one can generate new bound state solu- 
tions for various types of the well-known molecular potentials like the Rosen-Morse well 
[26] , the Eckart potential, the Hulthen potential [9] , the Woods-Saxon potential [5] and the 
Manning-Rosen potential [22] and others. It is also noted that under the PT-symmetry, the 
exponential potentials can be transformed into the trigonometric potentials with real bound 
state solutions. The KG equation with equally mixed scalar and vector Rosen-Morse-type 
potentials can be solved exactly for s-wave bound states. In the relativistic model, the 
energy equations for these potentials are complicated transcendental equations [26]. The 
non-relativistic limits are obtained with a proper replacements of parameters and/or by 
solving the original Schrodinger equation. The relativistic and non-relativistic results are 
identical with those ones obtained in the literature through various methods. 
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FIG. 1: Plot of the tRM potential [see Eq. (83)] for a set of parameters a = 0.5 and b = 17.0. 



TABLE I: The s-wave energy spectrum of the equally mixed scalar and vector Rosen-Morse-type 
potentials. 

n a qViV2 M El E2 E4 



1111- 


■1 4 


1.8137'' 


-1.9140° 


-3.3923 


" -3.9088 


2 




-2.2117 


-3.6791 


— 


— 


3 




-0.6606 


-3.3105 


— 


— 


4 




0.8879 


-2.7697 


— 


— 


5 




1.8766 


-1.9765 


— 


— 


1112- 


■2 5 


0.9989 


-3.7763 


-4.7275 


-4.9351 


2 




-4.1746 


-4.7795 


— 


— 


3 




-3.3814 


-4.5376 


— 


— 


4 




-2.3989 


-4.2008 


— 


— 


5 




-1.3083 


-3.7529 






1 0.5 1 1 


-1 4 1.9558 


-3.5288 


-3.8460 


-3.9773 


2 




1.9608 


-2.5367 


-3.5326 


-3.9216 


3 




1.2294 


-0.5126 


-3.0732 


-3.8358 


4 




-2.4823 


-3.7191 






5 




-1.7822 


-3.5695 






11 0.5 1 


-1 4 1.5783 


-3.2245 


-3.6502 


-3.9258 


2 




1.9995 


-1.5367 


-2.9520 


-3.7496 


3 




-1.9529 


-3.4736 






4 




-0.7335 


-3.0839 






5 




0.5489 


-2.5528 







■"Same as in Ref. [39]. 
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